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SUMMARY 

Vinti has recently modified his spheroidal potential so 
as to permit exact inclusion of the effects of the third zonal 
harmonic of the planet’s gravitational field. This cor- 
responds to a potential fitted exactly through the third zonal 
harmonic and about two-thirds of the fourth. The present 
paper treats the method for obtaining the position and ve- 
locity coordinates of a satellite moving in the field corre- 
sponding to this accurate modified potential. 
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COMPUTATIONAL PROCEDURE FOR VINTI’S 
ACCURATE REFERENCE ORBIT WITH INCLUSION 
OF THE THIRD ZONAL HARMONIC 

by 

N. L. Bonavito 


INTRODUCTION 

Vinti (Reference 1) has found a gravitational potential for an axially symmetric planet in oblate 
spheroidal coordinates. This solution accounts for all of the second zonal harmonic and more than 
half of the fourth zonal harmonic. This potential which simultaneously satisfies Laplace's equa- 
tion and separates the Hamilton Jacobi equation, succeeds in reducing the problem of satellite 
motion to quadratures. 

More, recently however, Vinti (Reference 3) has generalized his potential by means of a 
metric preserving transformation of the associated Cartesian system. This preserves separability 
of the problem of orbital motion when the potential coefficients J 2 a and J 3 are taken into account 
(Reference 4). The inclusion of J 3 is of considerable practical importance, permitting a more 
accurate treatment than that given by perturbation theory. This leads to the computing procedure 
for obtaining the position and velocity coordinates of a drag free satellite from a knowledge of its 
initial conditions. 


STATEMENT OF THE PROBLEM 

If we take r e as the earth's equatorial radius, and if, 

c 2 = r e 2 J 2 (i-? j 3 2 j; 3 ). 

S = -^ r eJ 2 ' 1 J3 - 

then the gravitational potential, 

V = - jU.(p 2 +c 2 17 2 )" 1 (p +778) 
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leads to a separability of the problem of satellite motion. Here, 8 % + 7 kilometers for the earth 
and the above potential leads to a fit of 

i 

I 


V 




Jn P n (Sintf) 


exactly through the third zonal harmonic, about two-thirds of the fourth zonal harmonic, and negligi- 
ble values of order j 2 3 for the higher harmonics. 

If we take p, v, and as oblate spheroidal coordinates satisfying the equations 


x + iy - r cos e 1!? 

= [(p 2 +c 2 )(l-r,a)] ,/ 2 e* . 

and 

z = r s i n £ = prj - 8 , 

then x, y, and z are the rectangular coordinates of a satellite in a Cartesian frame, with the origin 
at the center of mass of the planet. Also, r , 6 , and 4> are the planetocentric distance, latitude, and 
right ascencion. 

From Vinti (Reference 3), if a is the energy of the system, a 3 , the z component of angular 
momentum, and a v the separation constant, the generalized momenta are given by, = a 3 , 

P p = ± (p 2 +C 2 )- 1 F 1/2 (p) 

and 

Pt) = ± (l -p 2 )" 1 G ,/2 ( 17 ) . 

Here F(p) and G(i 7 ) are the quartics 

F(p) = c 2 a 2 + (p 2 + c 2 ) (- a 2 2 + 2pp + 2a j p 2 ) , 


and 


G(v) - ~ a 3 + (l _ V 2 ) { a 2 + 2m 7 ? 8 + 2o.j c 2 tj 2 ) . 
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The Hamilton- Jacobi function W(p, 77 , 4>) is then, 


W - + j*Ppdp + J*i 


P 77 d 7 ? 


or 


W = a 3 0 + 


f + (p 2 + c 2 ) 

J P. 


" 1 F 1/2 (p) dp + 


f± (l-^ 2 ) 

J n 


" 1 G 1 '' 2 (7 ? )d7 ? 


If /3j, ,fl 2 ' and ft ^ are constants of the motion, the orbit is then given by 


, <9W 

4 + ^ = 


0 . = -i*. 

^2 da. 


and 


, _ aw 

^3 da 


From Vinti's solution of these equations (Reference 3), together with the expressions for the gen- 
eralized momenta, we can describe a computational procedure similar to that described in Ref- 
erence 6. 


COMPUTATIONAL PROCEDURE 

Enter the initial conditions x. , y . , z. , x. , y . and z . for a time t . , with the constants p, r e , 
J 2 , J 3 < 0, and, 


J 3 2 J 2 " 3 ) s = -T r eJ 2 ' 1 J3 > 0 


(1.1) 


Compute, 


r. - Vx. 2 +v 2+ z. 2 i r.r. - x.x.+v.v.+z.z. 


( 1 . 2 ) 


r 2 + 2z. 8 + 8 2 - c 2 

p 2 = H < 1 + 


/ 4c 2 (z. + 8) 2 

V (r 2 + 2 z. 8 + S 2 _c 2 ) 2 


(1.3) 
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( 1 . 4 ) 


2 (*, + S) 2 

V - , (the sign of 77. - sign of (z. + 8) 

1 2 11 

Pt 


(the sign of 77. = sign of z ; ) 


Compute: 


where 


p. = _L (r . ^ + 8j.) + (T P' + + 2z i 8 + S2 ~ c2 > + 2c2 < z i + S > z 

2 Pi 


'|/( r ? + 2z. 8 + 8 2 ^ )~ 2 "+ 4 c 2 ( z.+ 8^) 2 


( 1 . 5 ) 


Vi = 


2c 2 t) 


(l'.T. + 8 z . ) + 

x 1 1 1 ' 


(r.r. + 8 z . ) (r 2 + 2z. 8 + 8 2 -c 2 ) + 2c 2 (z. + 8) z 
}^r 2 + 2 z i 8 + S 2 - c 2 ) 2 + 4 c 2 (z. + 8 ) 2 


( 1 . 6 ) 


-! = ^Ui 2 - pPi + ^i SH^i 2 ^ 27- ' 2 )' 1 


( 1 . 7 ) 


U 7 - x 2 + y 2 + z 2 

i i *7 1 i 


and 


Then, 


x i y ; - Vi x i 


(l -r?! 2 )- 1 [(p, 2 + c2 T 7 i 2 ) 2 V + a 3 2 - (l- 7 ?i 2 )( 2 a 1 c 2 7 7i 2 + 2 ^ 7 ); 8)] 


1 JL 

2 a, 


1 + 


2 a i a 2 2 


1/2 


J o( 1_e o 2 ) - and % 


( 1 . 8 ) 


Prime Constants 


Compute: 


X D 2 = - 2a 1 a 2 2 p- 2 and X D 4 


( 2 . 1 ) 
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P n 2 and Y n 2 = ly 


K o ~ Ti and K o 2 
p o 


(^1 + ^) = 2Po X d 2 [l-KoX D 2 Y D a -Ko 2 X D 2 Y D 2 (2X D 2 - 3X D 2 Y D 2 - 4 + 8Y D 2 )] 


( 2 . 2 ) 

(2.3) 

(2.4) 


P 0 2 X d 2 [i+K 0 Y d 2 (X d 2 -4)-K 0 2 Y d 2 (12 X d 2 -X d 4 -20X d 2 Y d 2 -16 + 32 Y d 2 +X d 4 Y d 2 )] (2.5) 


P 1 + ^2 


( 2 . 6 ) 


jgi ^2 


(2.7) 


/r 


( 2 . 8 ) 


-2 - 


a 2 2 - 2a j c 2 
2 ( a 2 2 " a 3 2 ) 


1 + 


1 + 


8 a i c 2 (a 2 2 


! )‘ 


(a 2 2 -2a lC 2 ) 2 j 


1/2 


§ - (sin 2 I ) - Tj Q 2 


(2.9) 

( 2 . 10 ) 


Mutual Constants 


p - a 


(l - e 2 ) 


(3.1) 


g a 2 c 2 f c 2 1 

" 2ac 2 (ap-c 2 §)(!-§)+ — p~ * 2 jl + — (3S-2)j§(l- 


S) 


(ap-c 2 ) (ap-c 2 s)+4a 2 c 2 S + — y 5 2 (3ap - 4a 2 - c 2 )§( 1 - §) 


(3.2) 


B - (2a) 1 ( ap - c 2 ) A + c 2 


„ , _ lj 

a o - a “ 2 A ’ p o 


(b + ap - 2Aa - c 2 ) 


~a~ r ~ 

a o 


2 = (Wo') 


1/2 


(3.3) 

(3.4) 


2_S_ X 2 
0 


Pn I (1_S) l 1 ' S"P7 S 


[l+ (c 2 /a 0 'p 0 ') ( 1 - 2S>] 


u" 1= i + (1717) (i-S) + ? 


(3.5) 
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(3.6) 


c 2 a I c 2 S A*2S W c 2 U \ 

^7p 7 u ’ C 1 ' V~^TpJ u ) Po' u V 1 " %'Po) 

/ c 2 SU\ 8 

p = l 1 " w)lV U(1 " S) 


= p + (p 2 + s ) 1/2 , 


P- (p 2 + s ) 1/2 


• ' - 




Using S', we now repeat steps (3.2) through (3.9) to obtain the quantities, 


A. B, a Q \ p Q ' a 2 ', e, U, C 2 , C r P, Vo - a" 0 * s 


(3.7) 

(3.8) 

(3.9) 


(3.10) 


Note: Equation (3.9) for S' together with the one step iteration (3.10) has now provided for us 

an accurate value of the element s that was originally approximated by Vq . (Equa- 
tion 2.10) 



B 1/2 . 



Using formulas of Reference 6, pages twelve to fifteen, we compute, 


= (1- 


e 2 ) 1/2 p 


X ^/ b _2 
p 


n = 2 


3 n(u;) R n-2 [(l -e 2 ) 1/2 ] 


= (1- 


e 2 ) 1/2 p 


_i £(^) p -G£H [( i - e2 ) i/2 i 

n=0 


(3.11) 


(3.12) 


where P n (b 1 /b 2 ) is the Legendre Polynomial of degreen, R n (x s ) = X"P n (x s 1 ) is a ploynomial of 
degree (n/2) in x s 2 , and x s = (1 - e 2 ) 1/2 . 

If m is an even integer compute, 


D 


= D 


2 > 
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If m is an odd integer compute, 


Then, 


D m = D 2i+1 


!>>-" Ip 


c /b, 


P 2n + llTT 


n=0 


co 

m= 0 


A n " 4 : (l - e 2 ) 1/2 p 3 e(- 2 b t b 2 2 p + b 2 4 ) 


A 12 “ 32 U - e2 ) 1/2 b 2 4e2 p“ 3 


2 1 


'22 


2 3 


(l - e 2 ) 1/2 p~ 1 e 


b x P 1 + (3b L 2 - b 2 2 )p 2 “ f bj b 2 2 ( 1 + ^7 ) p ~ 3 + |* b 2 4 (4 + 3e 2 ) p~ 


r ( 3b i 2 -b 2 2 )p _2 - | e 2 bjb 2 2 p- 3 + J 2 h 2 (6e 2 +e 4 )p- 4 


(l_e 2 )l/ 2 p -i 

(l_e 2 )l/ 2 p-i ^ (- bl b 2 2 p - 3 + b 2 4 p- 4 ) 


A 24 = 256 (l- e 2 ) 1 / 2 p-=b 2 4e4 


A 3i ~ (l-e 2 ) 1/2 p- 3 

A 32 = (l-e 2 ) 


2 + bj p~ 1 (3+ |e 2 )-p- 2 b 2 2 + c 2 ) (4 + 3^ 2 ) 


1/2 p- 3 


e 2 3 


1 „2 _„-2 I +c 2 )f4 e 2 + 


4 + 4 WP e J -P'^ \T 


2 e T 4 


'33 



CN 

J2 

CN 

1 

X 

(l_e 2 )i/ 2 -e 3 

- V u TJ 


A,. = - S + 

Q = (P 2 +S)' 


| 1/2 


( 3 . 13 ) 


7 



Then through third order in j 2 : 


= 1 - TC,p + (I c 2 + !c 2 )(p 2 + 4q 2 ) + M c 2 2 Q 4 


I Cj C 2 PQ 2 + -R5 C 2 C 2 Q 4 + ^ C/Q 6 


(3.14) 




= i Q 2 + P 2 -JC.PQ 2 +|p 2 Q 2 + £ (4C 2+ 3C 1 2 )Q 4 


” TF c i C 2 P( 5 4 + 235 ( 6C 2 + 15C, 2 C 2 )Q 6 + 2 048 C 2 3 Q 8 (3.15) 


- t c 2 " I c t 2 - (^c 2 c 2 + 4c 2 2 )(i + 4q 2 + 1 Q 4 ) - ^ C i C 2 


(3.16) 


1 - S’ h i " 2 


J (l +C, "C 2 ) _I/2 , h 2 = \ (1- C 1-C 2 )- 1/2 (3.17) 


e 2 = Q(l-P)" 1 , e 3 = Q( 1 + P) _ 1 


(3.18) 


277 V, = (- 2 a,')‘/* (a + bj+A^c'AjB^B - 1 )' 1 . 


277 v 2 = a 2 'U- 1/2 A 2 B 2 - 1 (a + b 1+ A,+c 2 A 2 B 1 'B 2 - 1 )- 1 


where e' < e < 1 


a + b. 


(3.19) 

(3.20) 


sgn a 3 [/ 2 p 0 ' (1 -S)] 


1/2 


I" T 7 


c 2 S \Po / \ a o Po ) 


1/2 


a n Pfi 


[ 


1 + aTPo' ( 1 " 2S) 


] 


(3.21) 


Here sgn a 3 £ o for a direct or retrograde orbit. 


Jacobi Constants 


B„ = 2PQ-8C,Q 3 


(4.1) 
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(4.2) 


B 12 = - (^ + 1 C 2 Q 4 ) 


B 13 

= Cj Q 3 / 24 

(4.3) 

B 1 4 

Q C 2 Q 4 / 64 

(4.4) 

B 2 1 

= c 2 p Q--^c lC2 Q 3 

(4.5) 

B 2 2 

= - 32 [(4C 2 +X 1 2 )Q 2 +3C 2 2 Q 4 ] 

(4.6) 


1 


B 23 

- -48 c i c 2 Q 

(4.7) 

B 2 4 

256 c 2 2 Q 4 

(4.8) 


COS <p. 


}[p} + C 2 y 1 - T?i 2 


sin0 i 


y^i 2 + c 2 / 1 - t ?. 2 


h 4 


p. 2 + T). 2 C 2 
p . 2 + c 2 


h 2 

V 


p. 2 + 7). 2 C 2 

1 — 77 i ' 2 


cos E . 


sin E 




^i h /.(^i 2 + c2 ) 

2a j' (p £ 2 + Apj + B j 


Vi - P 

sin ^i - Q 


cos i 


( 1 " ^i 2 ) 


/ ( a 2 2 “ a 3 2 ) _ 

7 S (l +C l^i" C 2 7 ?i 2 ) 


(4.9) 


(4.10) 


(4.11) 


(4.12) 
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cos v. 


cos E. - e 
1 - e cos E 


sin v . = 


(l - e 2 ) 17 2 sin E i 
1 - e cos E 


( 4 . 13 ) 


sin nv i for n - 2,3,4 
sin ro //. for n = 2, 4 


cos 30 


cos E^ 


e 2 + cos 


1 + e 2 cos 


♦ i) 
(*.*;) 


sin Ej'j 


(l - e 2 2 ) 1/2 sin (sij + -j) 
1 + e 2 cos («//. + "2 ) 


cos E 3 ' 


e 3 + cos 


1 + e 3 cos \\p. 


(*. - 5 ) 


sinE 3i 


(l-e 3 2 ) 1/2 sint^j - %) 

l + e 3 cos (> Ik ~\) 


2 'iUi + l) E 3i(^i-l) 


2 /Ttlf 


Xn 2 / 1 - 2£ 2 y 1 + 2& 

[-2a/)' 1/2 [b 1 E i +a(E i -esinE,) + AjV, + A n sinv. + A 12 sin 2 v.] 


+ a 2 


B' 0. + B n cos i//. + B 12 sin 2p i + B, 3 cos 30 £ + B 14 sin 4/n J t ; 
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/V = - a‘ (-2a/) 


- 1/2 


A 2 V 


4 

■ *E 

n= 1 


A,,,, sin nv. 

2n 1 


+ u 1/2 


(^i + t) B 2 +B 21 COS1 /'i 

B 22 sin 24 J i + B 23 cos 30 £ + B 24 sin 4^ i ] (4.18) 


/3 3 - 0. + c 2 a/ ("2a/)- 1 / 


1/2 


4 

A 3 v. + 2^ A 3] sin j Vi 

i=i 


u l/J {(l-S)-«[(h 1 .h 1 ) x , 1 


( h i h 2 )v u ] + B 3 ^ + ^CjCjQcos^. + ^ C/Q’sin^j (4.19) 


^ = ^'-u^b 2 | 


(4.20) 


c 2 /5,B/ 

= 2™i “ a/B 2 


(4.21) 


l o + go 


277 V . 


A 


Va/ /( a + b i + A i) A 2 1 


Orbit Generator 


(4.22) 


M s - l Q + 2n v 1 t 

= l o + Bo + 277 ^2 t (5.1) 


By Newton-Raphson iteration we solve M s +E 0 - e' sin (m s +E 0 ) = M s where M s +E 0 = 8; therefore, 


[8 n -e' Sin8 n -M s ] ! 

'8 n -e' sin 8 n - M s " 

2 

e ' sin 8 

n 

(l - e 1 cos 8 n ) 2 

1 - e ' cos 8 

n 


1 “ e ' cos 8 

u n 


and 8 n - M s initially. 


(5.2) 


cos v' - (cos 8 - e) ( 1 - e cos 8) 1 

sin v* = (l - e 2 ) 1/2 ( 1 - e cos 2)“ 1 sin 8 (5.3) 


11 





< 

O 

II 

< 

1 

s: 

(5.4) 



'p o = a 2 ' ( _ 2a / ) “ 1/2 U” 1/2 A 2 B/ 1 v Q 

(5.5) 

Mi = 

- (a+bj)' 1 

(Aj + c 2 A 2 B/ B/ 1 ) v q + ^ (-2a/) 1/2 U 1 ' a B 12 sin(2/-. + 2/' 0 ) 

(5.6) 



Mi e' M 1 2 sin(M. + E 0 ) 



E i " 

1 - e' cos (M + E n ^ 2 r , , x -i3 

^ s [l-e* cos(M s +E 0 )J 

(5.7) 


cos v" 

= j^cos (£ + Ej ) - ej |^1 - e cos ( e + Ej ] J 



sin v" 

= (l-e 2 ) 1/2 ^1 - e cos ( £ + Ej sin(£ + Ej) 

(5.8) 



v, - v" - (v Q +M S ) = (v* - v') 

(5.9) 

-B 22 B 2 - 1 sin(^ s + 2^ 0 ) + 

a 2 (-2a/)‘ 1/2 IT 1/2 B/ 1 [A 2 v 1 + A 21 sin (M s + v Q ) + A 22 sin (2M S + 2v 0 )] (5.10) 

- (a + b,)- 1 

A 1 v i + A u sin(M s +v 0 ) + A 12 sin (2M s + 2v 0 )+ (-2a/ ) " 2 U 1/2 {b/ 


+ B n 

COS^ s +^ 0 ) + 

2B 12 'P 1 cos (2 0 S + 2/< 0 ) +B 13 cos (3 4> s + 3 </< 0 ) +B, 4 sin (4^ s + 4>A 0 )j 

(5.11) 

_ 



m 2 




E 2 ' 1 - e ' cos (M s + E 0 + Ej) 

(5.12) 


II 1 _ 

cos v 

^cos (£ + Ej + E 2 ) - ejj^l - e cos (£ + E, + E 2 )j 



sin v' ,# = 

(l -e 2 )" 2 [l -ecos (£ + E, + E 2 )] _1 sin (£ + E, +E 2 ) 

(5.13) 



v 2 = v'" - (v 0 +M S + Vj) = (v'" -v") 

(5.14) 

-B 2 -> [b 21 cos (0 s +0 o ) + 2B 22 cos (2^ s + ^ 0 ) + B 23 cos (-*/>, + ^ 0 ) + B J4 sin (4^ s + 4"/> 0 )] 



+ a/ IT 1/2 (-2a/)' ,/2 B 2 - 1 [a 2 v 2 + A 21 Vj cos (M s + v Q ) + 2A 22 v, cos (2M s + 2v, 

i) 



+ A 23 sin (3M 3 + 3v 0 ) + A 24 sin (4M S + 4v Q 

)] (5.15) 
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Then, 


and 


E = 8 + Ej + e 2 , 


v “ M s + v 0 + V 1 + V 2 


'p - p s + p 0 + '/'l + p 2 

(*♦ 5 ) 


cos E 


e 2 + cos 


1 + e cos 


(** 5 ) 


(5.16) 


and. 


sin Ej' 


(l e 2 2 ) 1,/2 sin yp + 2 ) 


1 + e 2 cos 


(**t) 


cos E 3 ' 


e 3 + cos 


(♦-■») 


1 + e 3 cos - -j) 


and 


sinE 3 ' 


(1 -e/) 172 sin (</>-^) 
1 + e 3 cos ^ 

E 2 + "2 ) E 3 (^~ 2) 


x° 2 a- 21, 

2 yi + 24 

£ 2 - (0 + 3) 

E 3 (^- 5) 

Xi 2/1-24 

2 pi + 24 

p = a( 1 - e cos E) , 

T) - P + Q sin yp 


(5.17) 


(5.18) 
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and, 


<t> = Pi ~ ° 2 a 3 ( _ 2 a i') 


1/2 


4 ^ / 

A 3 v + £> SJ sin j y + 17 1,1/2 {(!" S)- 1/2 [( h i +h 2 ) 
j=l 


*0 


( h i “ h 2 )xi] +B 3^ + J C i C 2 Q cos '/' + -J2 C 2 2 Q 2sin ( 5 - 19 ) 


and 




p 2 + ^2 c 2 

/5 2 + C 2 


h 2 


V 


p 2 + TJ 2 C 2 
1 - T } 2 


h/ = (p 2 + c 2 )(l-p 2 ) 


(5.20) 


and 


ae ]/- 2a. j' (p 2 + Ap + B) 

P = 2 T~i 7\ sinE 

h P (p 2 +C 2 ) 


,2 -“ 3 ' 2 ) 


V 2 S ~~ ( 1+C 1 * " C 2 7)2 ) 

h 2 (l-T} 2 ) 


COS \jj 


)/ (p 2 + C 2 ) (l - T? 2 ) COS 0 


Y = y(p 2 + c 2 ) (l - 17 2 ) sin 0 


Z ~ PV 


X = XI 


PP 


p 2 +c 2 1-p 2 


Y- a 3 

V 


PP tit? 

Yl — — — I + 

p* + C Z 1 ~ Tj 2 


X ' a i 

V 


Z - PT) + PP 


( 5 . 21 ) 


(5.22) 
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REMARKS 


The accuracy of the orbit itself as a solution for the given potential is carried out through 
terms of the third order in j 2 , the coefficient of the second zonal harmonic. Its accuracy however, 
and thus that of the secular terms, may be increased at will. Periodic terms are carried through 
the second order, but their accuracy may also be increased. In order to carry this kind of accuracy 
perturbation methods presently in use would become far too cumbersome and impractical for orbit 
computational purposes. 

An advantage of accounting for J 3 in this way is the absence of small denominators in e or 
sin i that occur in perturbation theories. Thus, one can easily compute polar orbits and circular 
equatorial orbits. 

This program, similar in structure to the previous Vinti accurate intermediary orbit requires 
a relatively small number of storage locations throughout the entire computing procedure. Con- 
sequently, this program is expected to go as fast as the Vinti accurate intermediary orbit producing 
approximately 1800 minute vector points (time, x, y, z, x, y, and z, 1800 limes) each minute of IBM 
7094 Mod. I computer operation with simultaneous production of BCD tape output (Reference 6). 

Recent tests on the IBM 7094 have indicated that this latter Vinti method enjoys a rather 
sizeable advantage in computational speed over other methods which use perturbation techniques 
with the Vinti program producing some three to four times more minute vector points per unit 
time of computing machine operation. 

The residual fourth harmonic term (Reference 5) has been programmed and is presently being 
tested as is the J 2 j tesseral harmonic term. In addition, the orbit differential correction program 
written for the Vinti accurate intermediary orbit theory (Reference 7), will need only slight modifi- 
cation to account for the element S which replaces ! , (l r sin i in Vinti’s new theory. Consequently 
it is expected to go just as rapidly, producing a set of mean elements of even greater accuracy with 
which to predict satellite orbits over a longer interval of time. It has been shown (Reference 7) that 
the orbit differential correction program for the Vinti accurate intermediary orbit theory converges 
rapidly and with great accuracy. This latter program, tested on the Relay II Satellite with an ec- 
centricity of 0.23597617, and Satellite ANNA with an eccentricity of 0.00671710 gave the following 
results: 

Using radio direction cosine observation data, and for seventy equations of condition of the 
Relay II Satellite extending over a five hour arc following injection into orbit, the program converged 
on the third iteration to a standard deviation of fit criterion of 2.7 x 10" 3 within thirty seconds. 

Using Smithsonian Astrophysical Observatory optical data from the Anna Statellite extending over 
an arc of seventy five hours following injection, and with twenty equations of condition, the program 
converged on the second iteration to a standard deviation of fit of 0.2 x 10" 3 (milliradians). For 
seven equations of condition of Satellite ANNA covering the first forty-five hours following injection, 
the Vinti orbit differential correction program converged on the second iteration to a standard devia 
tion of fit criterion of 0.04 x 10~ 3 (milliradians). All tests were conducted using the IBM 7094 
Mod. I electronic digital computer. 
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